ON THE EXTINCTION PROFILE OF SOLUTIONS TO 
FAST-DIFFUSION 



PANAGIOTA DASKALOPOULOS* AND NATASA SESUM 

Abstract. We study the extinction behavior of solutions to the fast diffusion 
equation Ut = Au m on X (0, T), in the range of exponents m G (0, N ^ 2 ), 
N > 2. We show that if the initial data uq is trapped in between two Baren- 
blatt solutions vanishing at time T, then the vanishing behaviour of u at T is 
given by a Barenblatt solution. We also give an example showing that for such 
a behavior the bound from above by a Barenblatt solution B (vanishing at T) 
is crucial: we construct a class of solutions u with initial data uo = B (l+o(l)), 
near \x\ >> 1, which live longer than B and change behaviour at T. The be- 
havior of such solutions is governed by B(-,t) up to T, while for t > T the 
solutions become integrable and exhibit a different vanishing profile. For the 
Yamabe flow (m = jyTg ) the above means that these solutions u develop a 
singularity at time T, when the Barenblatt solution disappears, and at t > T 
they immediately smoothen up and exhibit the vanishing profile of a sphere. 

In the appendix we show how we remove the assumption on the bound on 
uq by a Barenblatt from below. 



1. Introduction 
We consider the Cauchy problem for the fast diffusion equation 



in the range of exponents < m < (N — 2)/N, in dimensions N > 3. The initial 
data Uo is assumed to be non-negative and locally integrable. 

Equation (jl.f [) arises as a model of various diffusion processes. It is found in 
plasma physics and in particular as the Okuda-Dawson low, when m = 1/2. Also, 
King studies (jl.ljl in a model of diffusion of impurities in silicon [8]. 
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When m = (N — 2)/(N + 2) equation Ijl.ljl is equivalent to the evolution of the 

4 

conformally flat metric gij — u N + 2 dxidxj by the Yamabe flow 

m = - Rg 

where R denotes the scalar curvature with respect to the metric g. The equivalence 

4 

easily follows from the observation that the conformal metric — u N + 2 dxidxj 
has scalar curvature 

.V - 2 

Au 7 ^ 
R = . 



u 

For an introduction to the Yamabe flow see in [TT] . 

Our goal in this paper is to study the vanishing behavior of solutions to the 
Cauchy problem (jl.l[) , under the assumption that the initial data uo satisfies the 
growth condition 

(1.2) u (x)=(j^J (l + o(l)), as\x\^oo. 

The self-similar Barenblatt solutions of equation (jl.ip given explicitly by 

C* (T-t) 



(1.3) B k {x,t) = 

with 



k(T-t) 2 t + |af 



N 
(1.4) 8= and 7 = - — 

satisfy the growth condition (|1.2p . The constant C* depends only on m and N and 
is given explicitly by 

(1 5) C * = 2m ( N - 2 - mN ) 

1 — m 

We will assume in the first part of this paper that the initial condition uo is 
trapped in between two Barenblatt solutions, i.e., 

. . / C*T \^ , . ( C*T 

(1 ' 6) {hTW) ^ W HvTN 

for some constants k± > lt2 > 0. As a direct consequence of the maximum principle 
we then have 

(1.7) B kl (x,t) <u(x,t) < B k2 (x,t) forO<i<T. 

In particular, u vanishes at time T. We will show in the first part of this paper 
that the vanishing profile of u is given by a Barenblatt solution. 
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Consider the rescaled function 
(1.8) u{x,t) = (T -t)- u(x(T -ty,t), r=-log(T-t) 



with (i and 7 given by (|1.4|) . It follows by direct computation, that u satisfies the 
equation 

(1.9) u T = Au m + H div(a? • «) 

and due to condition (|1.6|) . the inequality 



<L10 » (j^)^*^^ 

holds, for (a;, t) eR N x [- logT, 00). We denote by 

(1.11) B k{X ) = 

the rescaled Barenblatt solution. 

Our convergence results are described in the following two Theorems. The first 
result is concerned with the range of exponents < m < ~ 2 , for which 

the difference u — Bk, with Bk a Barenblatt solution, is integrable, amely we have 
J k jv {u — Bk) {x, t) dx < 00. Notice that this range of exponents includes the Yamabe 
flow, m=(N - 2)/{N + 2), when N < 6. 

Theorem 1.1. Let u solve the equation (|1.1|) for < m < N ^ 2 , with initial 
value uq satisfying 11. 6\) . for some constants k\,k,2- Then, the rescaled function u 
given by (|1.8|) converges, as t — > oo ; uniformly on R*, and also in L 1 (R ), to the 
rescaled Barenblatt solution Bk given by (jl.ll|) . for some ko > 0. T/ie constant ko 
is uniquely determined by the equality 

uo dx = / Bk dx. 

The second result deals with the range of exponents < m < -^5§, for which 
the difference u — Bk, with Bk a Barenblatt solution, is non-intcgrablc, namely 
J r n(u ~ B k ){x, t) dx = 00. 

Theorem 1.2. Let u solve the equation (II. ip . for < m < > N > 4, with 
initial value Uq satisfying H.b}) , for some constants k\,k2- Assume, in addition 
that 

(1.12) m = B ko +f 
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for some ko > 0, where Bk is a Barenblatt solution and f is in L 1 (M ). Then, 
the rescaled function u given by (|1.8p converges, as t — > oo, uniformly on K N , to 
the Barenblatt solution B^ . 

Remark 1.3. The condition < m < jrE^ m Theorem 11.21 implies that for any two 
Barenblatt solutions B^ and Bk>, B^ — B^ € L 1 (K Ar ) if and only if k = k' . 

One may ask whether condition (II. 2[) is necessary for Theorems 11.11 and 11.21 to 
hold true. We will show in section [5] that this is indeed the case. We will present 
an example of a class of initial conditions uo which satisfy the growth condition 

(1.13) u Q {x) = ( -j^p- J (l + o(l)), as \x\ oo 

with C* given by (|1.5p , for which the solution u of (|1.1|) with initial data ito satisfies 
the following Theorem. 

Theorem 1.4. There exists a class of solutions u of the Cauchy problem (|1.1[) wii/i 
initial data uq satisfying (| 1 . 1 3[) and wii/i f/ie following properties: 
(i) The vanishing time T* of u satisfies T* > T. 
(ii) The solution u satisfies as \x\ — > oo, the growth conditions 

2 

(1-14) u(x, t) > ) , on < t < T 

and 

(1.15) u(x,t) < C Kl , onT<t<T*. 

\x\ N - 2 

In particular, u becomes integrable on t > T . 
(Hi) The vanishing behavior of u is given by one of the self-similar solutions 
Q(x,t) (see section^ for the explanation of 9(x,t)). 

The vanishing behavior of the solution u in this case is described in the results 
of Galaktionov and Peletier [2], and del Pino and Saez [1] (in the case m = ]fe§), 
also formally shown by King [7]. 

The case of a special interest is m — , when the equation for u is eqivalent to 
the Yamabe flow of a corresponding conformally flat metric. The previous Theorem 
gives the following corrolary in the Yamabe case. 

Corollary 1.5. If u is a solution of the Cauchy problem (|1.1[) with initial data 
satisfying (JTTT3J) . Then the vanishing time of u is T* > T, (ii) in Theorem \ l-4\ 
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holds and 

(T* t)-r^u(x,t) - ( Cl )^, 
C 2 + \x - x\ 2 

as t -> T*, w/iere i£l lV and C x , C 2 > 0. 

Geometrically speaking, Yamabe flow starting at uq (described above) developes 
a singularity at T < T* at which the Barenblatt solution (cylinder) pinches off and 
immediately at t > T the solution becomes integrable. Due to the results of Del 
Pino and Saez it smoothens up at t > T and exhibits the behaviour of a compact 
sphere as t — > T*. 

The next section will be devoted on preliminary estimates for solutions u of 
with initial data in . The proof of Theorem 11.11 will be given in section [3] It 
will follow from the strong L^-contraction principle, Lemma 14. li which holds for 
the difference of the rescaled solutions u — Bf., for any Barenblatt solution Bf. . This 
method, based upon the ideas of Osher and Ralston [S], was previously used by 
S-Y Hsu in [5]. Since the difference u — Bf. ^ L 1 (K JV ) in the range of exponents 
< m < i^E§) f° r the proof of Theorem 11.21 we will need to weight the L 1 -norms 
with an appropriate power B^ . The proof in this other case is more involved and 
will be given in section 0J The last section will be devoted to the construction of 
the examples described in Theorem II .41 

Acknowledgements: The authors would like to thank S. Brendle, M. Gruneberg, 
J.R. King and J.L. Vazquez for helpful discussions. The second author would like 
to thank Max Planck in Golm (Potsdam) for the hospitality during the period in 
which part of this work has been done. 



2. Preliminary Estimates 

Our goal in this section is to establish the L 1 — contaction for solutions of 
and (|1.9p and some other preliminary results which can also be of independent 
interest. We begin by showing the following integrability lemma for the difference 
of any two solutions of (|1.1[) . 

Lemma 2.1. Assume that u,v are two solutions of (|1.1[) on l w x (0,T). If f = 
uq — vq G L 1 (R Ar ) and f is compactly supported, then u(-,t) — v(-,t) £ L 1 (K Ar ) for 
all t e [0,T). 
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Proof. Our proof is based on the well known technique of Herrero and Pierre [4] . 
We introduce the potential function 

(2.1) w(x,t)= [ \{u m -v m )(x,s)\ds 

Jo 

which satisfies the inequality 

(2.2) Aw > -|/|, on R N . 
Indeed, by Kato's inequality [6], we have 

A\ u m - v m \ > sign (u - v) A(u m - v m ) 
so from equation fTTTj) , we obtain 

(2.3) ^\u-v\<A\u m ~v m \- 

Integrating the previous inequality in time, and using that |/| = \uo — Vo\, we obtain 
(l2~2l . 
Let 

Z{X) = N(N-2). N L V^W^ lm ^ 
denote the Newtonian potential of |/|, so that from l|2.2[) we have 

A(w -Z)>0. 

Also, since |/| is integrable and compactly supported, there exists a constant C < oo 
for which 

(2.4) Z(x) < Vx » 1, 

where the range of x for which the estimate holds depends on the support of /. 
The mean value property implies that 

w{x) < Z(x) + — ^ [ (w(y) - Z{y)) dy 

UNP" JB p (x) 

< Z(x) -\ ^ / w(y) dy 

UnP JB p (x) 



for all x e R N ', p > 0. 
We next claim that 



(2.5) lim -4 / w{y)dy = Q. 
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Indeed, by Lemma 3.1 in [4] we have 

J ^ ^ \u-v\(y,t)dy < C M|uo - «o ||li(r«) + p N (J^j j 
which yields 

w(y,t)dx = \u- v\ m (y,s) dyds 

B p (x) Jo JB p (x) 



< C p 



t 

AT(l-m) 







(v m 
/ | m — w | dy \ ds 
Jb p (x) J 



(2.6) < C(T, 11/11) ^-2m/(l-m)_ 

Combining (|2.4[) . (|2.5[) . (|2.6p and letting p — > oo, we conclude the estimate 

(2.7) «;(!) < Z(a;) < j^^. 

We will use (|2.7[) along the lines of the proof of Theorem 2.3 in [4] to bound 
||(w — v)(-,t)\\ L imNy For < rjn < 1, let 77^ S Cfi°(18L N ) be a test function such 
that r?ij = 1 for \x\ < R and -q R = for |ar| > 2i?. Then |At?h| < C/R 2 and 
V?7i?| < C/R. Using equation (jl.ip . estimate (|2.7p and the integrability of /, we 
get 

\u-v\{;t)r) R dx < I \f\dx + 



f f \u m -v m \ A m dxds 

Jo Jb 2R \b r 



C f 

< II/IIlhr") + -m \ w(x,t)dx 

n Jb 2R \b r 

< C(\\f\\ LHMN) ) + C 

where C can be taken to be independent of R because of (|2.7I) . Letting R — > 00 in 
the previous estimate gives that 



sup / \u-v\(x,t)dx < C(||/|| £ i( H w)) +C < 00 
te[o,T) Jm N 

finishing the proof of the lemma. □ 

As a consequence of the previous Lemma, we will now establish the following L 1 
contraction principle for the solutions to (jl.l|) that are bounded from below by a 
Barenblatt solution B. 
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Corollary 2.2. Let u, v be two solutions of {TUj) with initial values uo, vq respec- 
tively and that f = uq — vq £ L 1 (M. N ). Assume in addition that u,v > B, for some 
Barenblatt solution B given by (II. 'S\i . Then, 

(2.8) / \u{;t)-v(-,t)\dx < / \u -v Q \dx, Vie[0,T). 

Jm N Jm n 

Proof. Let r/R be a cut off function as in the proof of Lemma l2.1[ with the support 

contained in B2R. Then, as before, we have 

^\ u (x,t)-v(x,t)\ < A\u m -v m \. 

If we multiply the above inequality by rjn and integrate over ]& N , since |A^i{| < j^, 
we get 

4/ \(u-v)(x,t)\r) R dx < / a(x,i) |(u - u)(M)l 
«e Jm N Jb 2R \b r 



with 



dO 



a(x ' t] = ] (e u + (i-9) v y-™ * cr2 > m B2R 



(since {9u + (1 - fl)?/) 1- ™ > B 1 -" 1 = C\(\x\ 2 + C 2 )). Hence, fixing t 6 [0,T), we 
obtain the estimate 

d_ 

dt 



\u — v\(x, t)r)n dx < C I \u — v\(x,t) dx. 

' B2r\Br 



\u — v\(x, t) dx < 

R N 



Assume that uq,Vq are compactly supported, so that / = uq — vq is compactly 

supported as well. Then, the right hand side of the above inequality converges to 

zero as R — > 00, due to Lemma |2~T1 This gives 

d 
dt 

which implies (|6 . 2[) . 

To remove the assumption that uq, vo are compactly supported, we use a stan- 
dard approximation argument. For any k > 1, we set Uq — «o XB fc (o)i v o — 
v oXB k (o)- Let u k , v k be the solutions of on H N x (0, 00) with initial val- 

ues Uq, Vq respectively. By standard arguments u k — > u and v k — > v uniformly on 
compact subsets of t™ x (0, 00). Also, by the previous argument 



\u"(-,t) - v k (-,t)\dx < / \uQ-VQ\dx< / \u -v \dx 
Jr n Jr n 

for all i > and for all k. Letting k — ► 00 we readily obtain (|6.2p . □ 

As an immediate corollary of Lemma |2. II we have the following result concerning 
the rescaled solutions u, v. 
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Corollary 2.3. Let u and v be as above. If uq — vq E L (R ), then for every 
t > — logT there is C(r) such that 

\u(x, r) — v(x, t)\ dx < C(r). 

Let u, v be two solutions of (jl.lj) satisfying (|1.6[) . We observe that the difference 
q = u — v satisfies the equation 

(2.9) q T = A(a(x, t) q) + | 7 | div(x • q{x, r)) 

on R N x [-logT, oo ), with 

Since both u and v satisfy (|1.10|) . for some constants ki,k%, it is clear that a(x,r) 
is smooth and satisfies the growth estimate 

(2.11) m(fc2 + N2) <a(^r)< m(fcl + N2) . 

Hence, (|2 .9[) is uniformly parabolic on any compact subset of x (0, oo). 
Let F(x,t) be a solution of 

(2.12) F t = A(ai{x,t)F) 
with 

i \ f 1 m 

a-\[X,t )— / ; ; ; dO . 

A direct computation shows that F(x, t) = F(x, (T — t)~< , t), r = — log(T — t), is a 
solution of equation 

F T = A{a(x, t) F) + | 7 | div(F • x) 
where a(x,r) is given by (|2.10[) . Similarly as before we have the following result. 

Corollary 2.4. Let F(x,t) be a solution of VTLJ\) . If F(x,0) e L 1 ^), then 
F(x,t) S L 1 (R N ). Moreover, F(x,t) e L 1 (R JV ) and for every r > - logT there is 
a C(t) such that 

\\F(;t)\\ L1{RN) <C(t). 
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3. The integrable case 

This section is devoted to the proof of Theorem 11.11 which deals with solutions 
of equation (jl.l[) in the range of exponents < m < In this case the 

difference of two solutions u, v satisfying (|1.7|) is integrable. We begin this section 
with the following strong contraction principle, which constitutes the main step in 
the proof of Theorem 11.11 Its proof as well as the rest of the argument is very 
similar to the proof of Theorem 2.3 in [5]. To facilitate future references we will 
sketch the proof of the strong contraction principle. 

Lemma 3.1. Let u, v be two solutions of (jl.ll) . for m £ (0, with initial values 

u ,v , satisfying [JTB\) . If 

mm{\\(u - w )+||l°°, II (wo - uo)-||z°°} > 0, 

then 

\\(u - u)(-,T)|| i i (R iv ) < \\u Q - v \\ L i(n."), t > -logT. 

Proof. Notice that by the comparison principle, u(x,t), v(x,t) satisfy (jl.lOp . The 
proof is almost the same as that of Lemma 2.1 by S-Y Hsu in [5], using the results 
in and [12]. 

Set q = u — v and observe, as above, that q satisfies equation (|2.9[) . Fix R > 0. 
By the standard parabolic theory, there exist solutions q^q 1 } of (|2.9[) in Qr = 
Br x (— logT, oo), with initial values <?(-, — logT) + , (?(•, — logT)_ and boundary 
values g+, g_ on 8Br x (0, oo), respectively. Notice that q+ — q 1 ^ is a solution of 
(|2.9[) in Qr, with initial value q(-, — logT) and boundary values g+ — g_. By the 
maximum principle, q = q+ — q 1 ^ on Qr. Similarly there are solutions q+ and q 1 ^ 
of (|2.9[) in Qr with initial values q + , and zero lateral boundary value. By the 
maximum principle, < q 1 ^ < q+ and < q 1 } < q^. Furthermore let qR be the 
solution of (|2.9p with initial value and lateral boundary value Bk 2 — B^ ■ By the 
maximum principle, we have 

0<q+,q*<q R - 

Let rj 6 C^°(R 2 ) be a cut-off function such that -q{x) = 1 on |a;| < 1/2, n{x) = 
for all \x\ > 1 and < rj < 1. Denote by vr = rj(x/R). The same computation as 
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in the proof of Lemma 2.1 in [5] gives 

\q(x,r)\r] R dx - \u -v\r] R dx 

( f a(x, t') q R (x, t') A Vr - | 7 | q R (x, t') x ■ Vf] R dx] dr' 

-logT \Jr n J 

a(x,T')q R (x,T') Ar) R - \-f\q R (x,T')x ■ V-q R dx ) dr' 



-logT \JR N / 

mm{q R (x, r'),q R (x, r')} r) R dx. 

The families of solutions q R (x,r) and q R (x,r) are monotone increasing in R and 
uniformly bounded above, which implies that 

q x = lim q R 

R — >oo 



and 



q 2 = lim q R 

R^oo 



exist and are both solutions of (1231) on R N x (—logT, oo ). This implies 
\q(x,T)\i] R dx - \u -v \r} R dx 



< ~5n I l a{x,T l )q R (x,T')dxdT l 

n -/-logT JR/2<\x\<R 



+ C / q R (x,T')dxdr' -2 / min{qf (x, r), q R (x, r)}ry H dx. 

J -\ogT J R/2<\x\<R JB Ro 

< C [ [ q R (x,T')dxdT' -2 [ mm{q R (x 1 T),q R (x,T)}r ]R dx. 

J — log T J R/2<\x\<R J B Ro 

By the same computation as in 5 , after letting R — > oo, we get 
\q(x,r)\dx— / | uq — vq \ dx 

R N JR N 



< -2 / mm{q 1 (x,T),q 2 (x,T)} dx 



3 Ho 

for all i? > and r > - logT. Since q x > qi Ro and q 2 > q 2 " , we obtain 



\q(x,r)\dx— / \uq — v\ dx 
Jr n 

<-2 mmlql^'ix^j^q^'ix^)} dx. 



Since qV* and q 2Ro are the solutions of (|2.9| in Q 2Ro with zero boundary value 
and initial values q+{-, — logT), q~(-, — logT), respectively, by the Green's function 
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representation for solutions, for any r > — logT, there exists a constant c(r) such 
that 

min q 2 + Ro > c(r) > and min q 2R ° > c(r) > 0. 

That finishes the proof of Lemma 14.11 □ 

The rest of the proof relies on the following result of Osher and Ralston ([H])- 

Lemma 3.2 (Lemma 1 in 9 ). Suppose that u(-,Ti) — > u in L 1 (R 7V ) ; as i — > oo, 
for some sequence t% — > oo. Let Bk be any stationary solution of ll.9\) . If v is the 
solution of \1.9\) in x [0, oo) with initial value v(x, 0) = u(x), then 

\\v(-,t) - -BfcHii^jv) = \\u- -Bfe||ii(RiV), Vr > 0, k > 0. 

Proof of Theorem We claim that in this integrable case there is a unique ko so 
that 

/ (m - Bk„) dx = 0. 

To prove that, let 

f(k)= f (u -B k )dx 

and observe that f(k) is a continuous, monotone increasing function with f{k{) > 
and /(fe) < due to (jT76j) . Therefore, by the intermediate value theorem there 
exists a unique ko such that /(fco) = 0. The rest of the proof is the same as in [5], 
based on the strong contraction Lemma 13. II 

4. The non-integrable case 

This section will be devoted to the proof of Theorem 1 1 . 21 which is concerned with 
the range of exponents < m < -^5f , N > 4. Assume that u is a solution of 
which satisfies the bound (|1.7p . Throughout this section u will denote the rescaled 
solution defined by (| 1 .8(1 , and Bk the rescaled Barenblatt solution given by . 

Since a difference of any two solutions u, v of equations is not always 

integrable in the range of exponents < m < t^E|> > 4, we need to depart in 
this section from the techniques used in [5] in which we heavily used the integrability 
of a difference of any two Barenblatt solutions. 

We define the weighted L^space with weight B a := (C*/(k 2 + |z| 2 )) Q /( 1-m ), as 

L 1 (B a ,R N ):={f\ [ f(x)B a (x)dx < oo}. 
Jr n 
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Lemma 4.1. Let u, v be two solutions of U.9\) . for < m < jrzti with initial values 
u 0l v satisfying ITb^) and u - v E L 1 ^", R N ), with a = (N — 2)(1 - m)/2 - 1 
(we will explain later such a choice of a). Let B := Bk 2 . If 

(4.1) max luo — «o| 7^ 

R N 

then, 

\\(u-v)(;t)B*ti r \\ l i < \\{u -io)B a 7 lR \\ L i, T > -logT. 
Proof. Set q — u — v. After rescaling, (|2.3[) becomes 

| ? | r <A(aM) + | 7 |V(a:-M), 

where a{x,r) is as in (|2~10|) . Let r) G C£°(IR 2 ) be a cut off function as before. 
Denote by t]r = r)(x/R) so that \Vtjr\ < C/R, \Ar] R \ < C/R 2 . Then, the above 
equation and integration by parts yield to 

\q(x,T)\rj R B a (x)dx - / \u - ii \r] R B a (x)dx = 
= f f (a{x,T , )\q\(x,T r ){B a {x)A m + AB a (x)r lR + 2\/ m WB a (x)} 

J- logT JR N V 

- \~f\\q\{xy)x- {B a {x)V m + VB a {x)r lR }dx}dT l . 
Moreover, 

\q(x,r)\r] R B a (x)dx - / \u - v \ rj R B a (x) dx 



<^S I I a(x : r')\q\(x,T r )B a (x)dxdT' 

- logT JR/2<\x\<R 



R 2 

c 



(4 2) ° J - logT JR/2<\x\<R 



a(x,T')\q\(x,r') |V B a (x)\ dxdr' 



+ C I / \q\ (x,T')B a (x) dxdr' 

- logT J R/2<\x\<R 



f I {a{x,T')AB a -\j\x-VB a }\q\(x,T')r) R dxdr' 

J- log T JM. N 

: h(R) + I 2 (R) + h(R) + h{R) V R > R > 0. 



We fix 



(A *\ (N~2)(l~m) 
(4.3) a = 1. 

Claim i. If max R « \u(x,t) — v(x,t)\ ^ 0, there exists a constant C(r) > 0, such 
that h(R) < -C(t). 



14 
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Proof of Claim i. We recall that 

m(|a;| 2 + fc 2 ) m (\x\ 2 + fc x ) 
^ < a{x, t) < — 

with C* = 2m ^\_?~ mN ^ ■ A direct computation shows that 

ABa _ 2 a ( [2a - (1 - m) (JV - 2)] \x\ 2 -k 2 (l- m) N) ^ < Q 
(1 - to) 2 (1 + |x| 2 ) 2 

provided that a < (1 - m)(N - 2)/2. Hence 

a(x, t') AB a (x) - | 7 | x ■ VB a (x) < m (|x J^ + k ^ AB a {x) -\y\x- VB a {x). 
Again, by direct computation (using also that I7I = jy-2-jVro anc ^ ^hat ^* ~ 

2m(N—2—Nm)\ n 1 

— —, '- we hnd 

1— m 

m {\x\ 2 + hz) 



C* 



-AB a (x) - \ 1 \x-VB a {x) 

k 2 (N-A-m(N-2))N 0(m,n,k 2 ) 

t> = 57 ; < U 



2{N(l-m)-2){k 2 + \x\>) {k 2 + \xf 



for m < ^— I and a = — 2 If- 1 — 1. From this the claim easily follows. 



iV— 2 2 

We will now compare the terms Ii in (|3.ip in order to get a strong contraction 
principle with the weight B a . 

Claim ii. There is a uniform constant C (independent of R) such that 
(4.4) / \q(x,T)\r )R B a (x)dx< [ \u - v \r) R B a (x) dx + C. 

JR« JR N 

In particular, if u - v E L 1 (B a ,'R N ), then u(x,t) -v(x,t) 6 L^B",^). 
Proof of Claim ii. The proof of Claim ii is similar to the proof of Lemma l2.ll once 
we know 



h{R) < C f f \q\(x 1 T , )B a {x)dxdT 

J- loeT J R/2<\x\<R 



logT JR/2<\x\<R 

r r dx 

< C / 2 <C(r), »€ {1,2,3}, 

J- logT JR/2<\x\<R \ x \ 1 ~ m 

for our choice of a, where C is independent of R (for l2(R) see also the arguments 
in the proof of Claim iii) . Estimate (|4.4[) now follows from the computation in (|4.2|) 
and the claim i. 

Claim iii. We have, limfl^oo Ii(R) = 0, i = 1, 2, 3. 
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Proof of Claim Hi. By the Claim ii we have |g|(-,r') e L 1 (B a ,R N ), for r' £ 
[— log T, r] . Hence, by choosing R sufficiently big, the integral 

f f \q\{x,T r )B a {x)dxdT' 

J-\ogT JR/2<\x\<R 

can be made arbitrarily small. This readily implies that Hdir-hx, J3 (i?) = 0. In 
addition, since a(x,r) < C\x\ 2 , this also implies that lim^j^oo Ii(R) — 0. Fi- 
nally, observing that |V_B Q | « B a / R on R/2 < \x\ < R, we conclude that 
liniR^oo Ii{R) — 0, finishing the proof of the claim. 

Combining the above two claims and (|4.2j) , concludes the proof of Lemma l4~Tl □ 

Since we have Lemma |4~T| the following result holds due to Osher and Ralston if 
we replace the usual L 1 norm by the weighted L 1 (B a ,M. N ) norm. This replacement 
will leave the proof of the following Lemma unchanged. 

Lemma 4.2 (Lemma 1 in [9] by Osher and Ralston). Let R > Rq. Suppose 
\\u(-,Ti) — vo\\ L1 (gc, RJV -j — > as i — > oo ; and let Bk be any stationary solution of 
S1.9\) . If v is a solution of U.9\) in M. N x [0, 00) with initial value v(x,0) = vq(x), 
then 

\\v(-,t) - B k \\ L1{S , aMN) = \\v - B k \\ L1( g a RN) , 
for all t > and all k > 0. 

Remark on the proof of Lemma \4.2\ Define T(t)uo — u(t), where u(t) is a solution 
of (If ,9p starting at uq. The proof of Lemma l4~2l uses only that T(t) is a semi-group 
on an L 1 (B a rjR, R Af )-closed subset of L°°, satisfying the contraction principle (|6.2|) 
and fixing Bk- 

The following simple convergence result will be used in the proof of Theorem 

n~2i 

Lemma 4.3. Let uq satisfy ll.6\) for some constants k%, hi- Take any Ti — > 00 and 
let Ui(-,r) = u(-,Tj +t). Then, passing to a subsequence, Ui converges, as i — > 00, 
uniformly on compact subsets of M. N x (—00,00) to v(x,t), an eternal solution of 
JO|) . satisfying (TW\) . 

Proof. Since u satisfies (|f .1 01) . equation (|f .9\i is uniformly parabolic on B2R x 
[— (logT)/2 — Tj,oo), for any R > 0. By standard parabolic estimates, the se- 
quence {ui} is equicontinuous on compact subsets of x (—00,00). Hence, by 
Arzela-Ascoli theorem and the diagonalization argument any sequence {Hi} will 
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have a convergent subsequence, converging uniformly on compact subsets to v, an 
eternal solution of (|1.9|) satisfying (|1.10|l . □ 

Recall that we denote Bk simply by B. 

Claim 4.4. The sequence Ui(x,r) converges to v(x,t) in L^B a ,R N ) norm. 

Proof. By Lemma 14.11 and Lemma 13.11 we have 

/ \ui(x, t) — B(x) \ dx < C, and 

\v.i(x,T) -B(x)\ ■ \B a (x)dx < C, 
where C is a constant independent of i and r and therefore, 

\ui(x,r) — v(x, r)\B a (x) dx < 
< I \v,i(x,T) — v{x : r)\B a {x) dx + 



B' 



+ (/ \u t (x,T) - B{x)\B a (x)dx+ / \v(x,t) - B{x)\B a (x) dx) 

J\x\>R J\x\>R 

< / \ui(x,r) - v(x,r)\B a (x)dx + 
Jb r 

+ mdni-m) ( / \u l {x,T)~B(x)\dx+ f \v(x,r)-B(x)\B a (x)dx)< 

n J\x\>R J\x\>R 

/C 
^ \ui(x,T) - v( Xl T)\B a (x)dx + R2a/(1 _ m) ■ 

If we let i — > oo in the previous estimate, since r) — » r) uniformly on 
compact sets, we get 

C 



lim / \iii(x, t) — v{x, r)\B a {x) dx < 



jf2a/(i-m) ' 

which holds for every R > and therefore ||uj(-,t) — v(-, T )\\Lim a (x) r n ) ~* as 
i —> oo. □ 

Proof of Theorem \1.2{ 

Proof. For a sequence Tj — > oo, let a subsequence for which u(-,Tj) — > {Jo, as 
ife — » oo, uniformly on compact sets of K , as shown in Lemma 14.31 We will show 
that vq = B, as stated in the Theorem. By the previous claim and Lemma 14.21 we 
have that 

(4-5) \\v(;t) ~ B(-\\ Ll{§a RN) = \\v(;t ) - B(-)\\ Ll{§a RN) , 
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for all r > —logT. On the other hand, if max \v(x, 0) — B\ > we have the strong 
contraction principle (|4.1j) . which contradicts (|4.5|) . Therefore, v(x f 0) = S(a;). □ 

5. Solutions that live longer 

In the previous sections we established the vanishing profile of solutions u of 
equation with inital data uq trapped in between two Barenblatt solutions 

with the same vanishing time T, i.e. when (|1.2|) holds true. We showed that if 
Mo satisfies (|1.6p . then u vanishes at time T and the rescaled solution u(x, r) = 
(T — t)~@u(x (T — t) 7 ,t), with r = 1/(T— t), converges, as r — > oo, to a rescaled 
Barenblatt solution B. 

In this section we will show the condition (|1.6| is necessary for Theorems 11.11 
and 11.21 to hold true. We will prove Theorem 11.41 which presents an example of a 
class of initial conditions uq which satisfy the growth condition 

( C* T\ 

(5.1) u (x) = T12- + as|x|^oo 



with C* given by (jl .5|) , for which the solution u of (|1.1|) with initial data uq vanishes 
at time T* > T. In addition, we will show that the solution u remains strictly 
positive for t < T* and it satisfies, as \x\ — > oo, the growth conditions u(x, t) « 
C(t) \x\~t^, with C{t) > on < t < T and u(x,t) = 0{\x\~^), on T < t < 
T* . In particular, u becomes integrable on t > T and its vanishing behaviour is 
given by a class of self-similar solutions 9(x, t). 

It is well known that the Barenblatt solutions given by (| 1 . 3[) are not the only 
self-similar solutions of equation (II. ip . It was shown in [10] that (jl.ip possesses 
self-similar solutions of the form 

r 1 —10 

e(r,t) = (T-trf( V ), r, = -—- e ,a=- , 

(1 — ty 1 — m 

where the function / is a solution of an elliptic non-linear eigenvalue problem with 
eigenvalue 9 satisfying the bound 

m 1 

< 6 < 



(l-m)N-2 2' 

/'(0) = and f(rj) = 0(r]~^ N ^ 2 ^' m ) as 77 — » 00. The solution / was shown to be 

unique apart from a scaling due to the invariance of / under the transformation 

2 

f(r); A) = A~ 1 - m / (77/ A; 1) and satisfies 



N-2 



(5.2) /(»?;!) as 77 
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It was shown in [3j that for any radially symmetric solution of (11. ip with initial data 
satisfying the growth condition uo(r) = 0(r ~) as r — > oo, then the vanishing 
behavior of u is described by the self-similar solutions 0, i.e. there exists A > 
such that the rescaled solution satisfies 

(T-ty a u( V (T-t) B ,t)^ f( m \) 

uniformly in rj > 0. 

In the proof of Theorem 11.41 we will use the following lemma, which is also of 
independent interest. 

Lemma 5.1. Assume that v is a solution of (|1.1[) on M. N x (0,T). Assume that 
t'o < / with f G L 1 (]R Ar ) and radially symmetric. Then, at time t > the solution 
v satisfies the bound 

C 

v(x,t) < | |JV _ 2 , \x\ > 1. 

Proof. We introduce the potential function 

w(x, t) = v m (x, s) ds 
Jo 

which satisfies the inequality 

(5.3) Aw > -/, on R N . 
Let 

Z[X) = N{N-2)u N L V^W^ HV) dV 
denote the Newtonian potential of /, so that we have 

A(w -Z)>0. 

Since / is integrable and radially symmetric, there exists a constant C < oo for 
which 

(5.4) Z{x) < Vx 6 



<>N 



Indeed, this follows from the observation that for a radially symetric / the Newto- 
nian potential of / is also given by 



ZU-). I -±zt I f(y)dyd P . 
P J\y\<P 



Similarly as in the proof of Lemma 12.11 we get 

C 



(5.5) w(x)= / v m (x,s)ds < Z(x) < 
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We will now use this bound together with the Aronson-Benilan inequality 

(5-6) v t < T 1 v 

(1 — m) t 

to conclude the desired bound on v. Indeed, we first integrate (|5.6p in time to 
obtain the inequality 

v(x,t 2 ) . ft 2 \ jA ^ ... . 4 
7 < r ' if < <i < t 2 . 

V{X,ti) \tij 

Hence 



w(x)> v m (x,s)ds>v m (x,t) (|V "* = c(t)v m (x,t) 

.It lo .It lo^tJ 



It/2 Jt/2 

which combined with (|5.5p implies the bound 

C(t) 



v(x,t) < 



N-2 
X\ m 



as desired. □ 

We now proceed with the proof of the Theorem. 
Proof of Theorem \l-4\ We begin with the following observation: If w is a solution 
of (jl.ip which vanishes at time T, then for any a > 0, the solution of given 
by W(x,t) = w(ax,a 2 t) has vanishing time T 1 = T/a 2 . Hence, we can make T 1 
arbitrarily large by choosing a sufficiently small. 

Let / > 0, be any radially symmetric integrable function such that f(x) 
. as \x\ ~ > oo. In particular, we can take f(x) to satisfy f(x) — 

N — 2 

0(\x\ ~), \x\ — » oo. Choose a sufficiently large so that the vanishing time T 1 of 
the solution w of (jl.ip with initial data wq(x) = f(ax) satisfies T 1 > T, with T as 
in fl5jj. Set 

1 1 

and let B(x, t) — (^jJ^y^J 1 ™ the Barenblatt solution with B(x, 0) = ^ i ^+i ^ 1 
Denoting by u the solution of (jl.ip with initial data no it is clear application of 
the comparison principle that u > w so that the vanishing time T* of u satisfies 
T* >T > T. This proves (i). Also, since f > 0, u > B, so that (flTil)) is satisfied 
as well. 

Since u - B(-, 0) = / g L 1 ^), by Corollary! 



[[«(-,*) -fl(-,i)|| il(R *) < ||/(-)||ii(R*). 



2(1 
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Since B(-, t) = for t > T, this implies 

(5.7) \\u(;t)\\ L1{m) < ||/|| L1(R «) for t G [T, T*). 

Combining the estimate (|5.7p and Lemma |5. 11 if we take v(x,0) = u(x,T), yields 
to (fl~T5l) . 

The statement (iii) of our Theorem now immediately follows by the result of 
Galaktionov and Peletier in [SJ. □ 
The proof of Corollary 11.51 now easily follows. 



Proof of Corollary \1.5\ It is known that if m — then a = ^±1 arLC j 6 = 

in the definition of self-similar solutions 0(ar, t). In this case function / is given 
explicitly by 



/(»?,A) 



N + 2 

K N X ' 



X 2 + ?y 2 / 

The proof of Theorem 1 1 . 41 implies that for t > T, we have the bound 

By the result of Del Pino and Saez in [T], if the vanishing time of u is T* with 
T* > T, then there is A > so that 

(T* - i)~T^r u (a; 5 1) _> /(|a;|, A) as i -> T*. 

□ 

We will end the paper with the following remark regarding the examples of 
solutions u constructed in Theorem 1 1.41 We know that these solutions u in live up 

N — 2 

to T* > T, and satisfy u(x,t) > B(x,t) for t < T and u{x,t) < C{t) \x\ — ~ for 
t > T. 

Remark 5.2. In the examples of solutions u constructed in Theorem \l-4\ the 
Barenblatt B(x, t) dies off exactly at time T , in the sense that there is no sequence 
\x%\ — > oo and A > so that 

(5.8) u(x h T)> A 2 . 



Proof. Lets argue by contradiction. Assume there is a sequence \xi\ — > oo so that 
pTB|l holds for some A > 0. Take t > T/2 so that 2 m (C*(T < A/2 and let 

5 = (T - i)/2. Notice that by our choice of /, < u(x,0) - B{x,0) G L 1 ^). By 
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the comparison principle, we have u(x, t) > B(x,t) for t > 0, and by Lemma 1 2. l\ 
u(x, t) — B(x, t) £ L 1 (R W ). Let 



(u m ~B m )ds. 



As before, Aw > -(u - B) and since \\{u - B)(-,t)\\ L(RN) < ||/||ii( R JV), using 
Newtonian potentials we obtain that 

for a uniform constant C. This yields the existence of s 6 (t, t + S) such that 

u™(x, S )<B™(x,s) + —^ 

(5.9) t ^ m 

/ C*(T-s)\~ C 



< 



5- \x\ N - 2 ' 

For each Xi choose Si E [t, t + S] so that (|5.9[) is satisfied. The Aronson-Bcnilan 
inequality gives 

T 

u(xi,T) < — u(xi, Si) < 2 u(xi, St). 

Si 

Combining (|5T5|) and ([53]) yields 



\Xi 

and therefore 



A C 

< 



which can not be fulfilled for \xi\ >> f . This finishes the proof of our claim. □ 



6. Appendix 

In this appendix we improve Theorem 1 1.21 The goal is to remove the assumption 
u (x) > B kl {x,0), where B kl (x,t) = (C* (T-t)/(ki + \x\ 2 ))^ is a Barenblatt 
solution with the same vanishing time as u. This was assumed in (|I.6[) . The 
bound from above is necessary, as was proven in the previous section. Denoting by 
B(x, t) = (C* (T - t)/(k + \x\ 2 ))^, for some k > 0, we will show the following 
result. 

Theorem 6.1. I/O < m < for N > 4 and 

(6.1) \u (x)-B(x,0)\<f(\x\)eL 1 (R N ) and u (x) < B(x, 0) 
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for a positive, radial function g, then u vanishes at the same time as B(x,t) and 
the reseated solution u(x,t) converges, as r — > oo, uniformly M. N , to the reseated 
Barenblatt B. 

To simplify the notation we will assume that B(x,t) = (C* (T - t) / (l+\x\ 2 )) T ^ . 
The proof of the theorem will be based on a sequence of observations. 

We begin by noticing that condition (|6.1j) implies the ^-contraction principle 

(6.2) / \u(;t)-B(-,t)\dx< I \u o -B(;0)\dx, Vie[0,T). 

The proof is the same as the proof of Corollary 12. 21 which goes through under the 
weaker assumption that only vq > B (and not necessarily uq > B). Furthermore, 
this implies that u and B have the same vanishing time. 

In order to be able to take the limit of the rescaled solution, we need to establish 
the necessary a'priori estimates. It turns out that it is possible to do so, just by 
using the fact that the difference u(-, 0) — B(-, 0) G i 1 (R Ar ) and that B and u have 
the same vanishing time. Introducing as in the previous sections the rescaling 

(6.3) u{x,t) = (T -t)-Pu(x(T -ty,r), r = -ln(T-t) 
(with /3, 7 given by (|1.4p ~), which satisfies the rescaled equation 

(6.4) (u) T = Au m + | 7 | div(a; • u), on R N x (- logT, oo) 
we have: 

Proposition 6.2. If \u (x) - B(x,0)\ < f(\x\) G L 1 (M Ar ), then there are positive 
constants Ci,C2,ro,To such that 

C C 

( 6 - 5 ) (r 2 + 1) l 1 /(l- m ) ^ T) ^ P + l)V(l-m) ^ " ^ r ° 

for all t G [to, oo). 

Proof. We will first prove (|6.5p under the assumption that uq is radially symmetric. 
At the end of the proof we will remove this assumption. 

As we have observed in the proof of Lemma |2~T| absolute value \u — B\ satisfies 
the differential inequality 

— \u-B\ < A\u m -B m \. 
at 

Hence, for any fixed numbers T/2<s<t<T, the function 

w[x) = / \u m -B m \(x,l)dl 
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satisfies 

Aw > -\u - B\(s). 
Let Z be such that AZ = — \u— B\(s) and therefore 

A(w - Z) > 

which, by the mean value property, as in the proof of Lemma 12.11 yields to the 
inequality 

w{x) < Z(x, s), r = 
Since, u — B is radially symmetric, the function potential Z is given by 



Z(x, s) 



I 17=1 / \u - B\(s) dy dp, 

Jr P J\v\<P 



\v\<p 

for an appropriate constant c n . We conclude that 



/' 



which combined with (|6.2|) and our assumption (|6.ip implies the bounds 

B m (x,l)dl-cM^ < f u m (x,l)dl< I B m (x,l)dl + C^^. 

^ J s J s ^ 

Next, fix t e [3T/4, T) and choose s S [T/2, T) such that T - t = t - s, so that 
r^(f-s)r(x,s)< / B m {x,l)dl <2^ (t- s)B rn {x,t). 



Il 1 



Combining the two last inequalities gives 



ll/IU* / i 



2"— B m (x,s) - C^^r < / u m (x,l)dl < 2— B m (x,t) + C- 

r" z t — s J s /- 

^From the above we conclude a pointwise bound from above and below for u m (x, t) 
with the aid of the Aronson-Benilan inequality, (log it)t < 1/((1 — m) t), which after 
integration implies the bound 

d)T^u m Cx,t) < — !— [ u m {x,l)dl < (-)t^ u m (x,s). 
t t — s J s s 

Due to our choice for s (for a given t) we have t/s < 2. Hence, combining the last 

two inequalities gives 

(6.6) u m {r,t) <C 1 B m {r,t)+ ' 



(T-t) r N ~ 2 

which holds for all t e [3T/4,T), since t is arbitrary, and also 

C 



(6.7) u m (r,s) >C 2 B m (r,s)- 



(T - s) r^- 2 
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which also holds for all s G [3T/4, T) since t is arbitrary and T — t = t — s. Rescaling 
inequalities ()6.6p and ()6.7|) . and using (|1.4p we conclude 



< u m (r,r) < 



(j.2 _|_ 2)m/(l-m) r W-2 — V'';— / j _|_ ]\m/(l-rn) r iV-2 

for some uniform constants C%,C2,C and t > tq := — log(T/4). This readily 
implies (|6.5[) for r >r$ (independent of r) if iV — 2 > j- 2 ^;- The last is equivalent 
to m < an( j is implied by our assumption m < (N — 4)/(iV — 2). 

In the case where u(x, 0) is nonradial and B(x, 0) — u(x, 0) is bounded from above 
by a radial function in i 1 (M Ar ), define 

u (r) := inf u(x,0) and Tio(r) := sup u[x, 0). 

By our assumption, B(x,0) — f(r) < u(x,0) < B(x,0) + /(r), where r = |x| and 
/ G L 1 (R 7V ). Since B(x, 0) is a radial function itself, we have B(x,0) — u Q (r) G 
L 1 (R JV ) and B(x,0) - u (r) G L 1 (M iv ). Let w(x,t) and u(x,t) be the solutions to 
(II. ip with initial data w ( r ) anc ^ "o( r )j respectively. Then, the radial result implies 
the bounds 

C\ C , . C2 C 

< u (r, t) < 



f r 2 _|_ ]\m/(l-m) r 7V-2 ' / — ^ r 2 _|_ ^m/(l-m) r W-2 ' 

and 

< u (r t) < 

W 1 iW(l-m) r JV-2 — V'!';— u 1 r N-2 " 

Since, by the comparison principle, w(r, r) < u(x,r) < u(r, r), the above inequali- 
ties imply (|6.5p in the nonradial case. □ 



Proof of Theorem \6.1\ The previous Proposition yields the existence of ro > and 
to < 00 so that 

Ci . C 2 

< u(r, r) < 



(r 2 + - v ' ' - ( r 2 + ^l/Cl-m) 

for t G [to, 00) and r > ro. Let Q ro = B ra x [ro, cx)). Hence, there exists a constant 
Co = c(ro, ro) > such that 

C 

" (r ' T) " (rg + l)"/(i-"') > 011 Wo ■ 
By the maximum principle 

Co 



inf m > — „ 
<2"o " (rg + l)V(i-m) 



which, combined with the lower bound in (16. 5p implies that 

C 

T ) ^ ( r 2 + 1 )l/(l-,n) ' 011 X t T 0' °°) 
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for a constant C\ that depends on ro- 

By our assumption we have u(x,t) < B(x) on M. N x [— logT, oo). Hence, there 
are uniform constants C\ and Ci such that 

( 6 - 8 ) p + iy/(i- m ) ^ fi fo r ) ^ ( r 2 + ■ onR^ x [r ,oo). 

We conclude that the difference u — B satisfies the equation 

(u - B) T = A(a(u - B)) + |7| div(a; • (u - £?)), 
with a(x, r) = /„ ^-^-jl^-^ satisfying the bounds 
(6.9) Ci (r 2 + 1) < a(x, r) < C 2 (r 2 + 1) 

on R N x [to, 00). The rest of the proof of Theorem 16.11 is the same as that of 
Theorem O □ 
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